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Abstract Acoustic waves propagation of in composite of 
water with embedded double-layered silicone resin/silver rods 
is considered. Approximate values of effective dynamical con- 
stitutive parameters are obtained. Frequency ranges of simul- 
taneous negative constitutive parameters are found. Local- 
ized surface states on the interface between metamaterial and 
"normal" material are found. Doppler effect in metamaterial 
is considered. Presence of anomalous modes are shown. 
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1 Introduction 

Media with effective negative constitutive parameters attract 
attention for several decades due to unique properties of wave 
propagation in them. In pioneering works of Veselago fTJ and 
Pendry [2| electromagnetic waves in media with simultane- 
ous negative e and p were considered. Since that time, a lot 
of theoretical and experimental works has been done on this 
subject [3]. In fact, for acoustic waves metamaterials are pos- 
sible too and for them both density and stiffness should be 
negative. Some proposals and experiments to achieve nega- 
tive effective constitutive parameters for acoustic composites 
were made in the past, exploiting different techniques, such 
as: embedding soft inclusions in fluids (4|5), using Helmholtz 
resonators [6] and pipe -membrane structures (7). 

The aim of this work is to show that it is possible to con- 
struct a metamaterial to exhibit negative constitutive param- 
eters within ultrasonic range and demonstrate that acoustic 
waves can exhibit exotic behavior, namely presence of local- 
ized states at the boundary between metamaterial and "nor- 
mal" material and reversed Doppler effect. In the mathemat- 
ical model both longitudinal and shear modes are included 
in calculations. Also internal dissipation losses are taken in 
account. 



2 Mathematical model 

2.1 Effective constitutive parameters 

In this work a composite consisting of fluid host with rods 
embedded in it is considered. Rods have round cylindrical 
shape and consist of elastic materials. Acoustic waves propa- 
gate in host material (density po and bulk modulus Bo) with 
wave-vector perpendicular to generatixes of the rods. A meta- 
material should be quasi-isotropic, so the wavelength is much 
longer than the distance between the centers of cylinders (L) 
and their radii (R) [3|. All cylinders are considered identical 



and are placed at random at approximate equal distance be- 
tween neighbors. The aim is to compute the dynamical effec- 
tive constitutive parameters of a metamaterial on given fre- 
quency and dispersion of the propagating bulk acoustic wave. 

To proceed the calculations coherent potential approxi- 
mation (CPA) [5.8] is used. We used this method applied 
to fluid/solid composite taking into account longitudinal and 
shear polarizations in solid inclusions and acoustic dumping 
in both inclusions and host. 

Following the CPA procedure, one inclusion is surrounded 
with cylinder of radius L (forming coated inclusion) and out- 
side the cylinder some effective fluid is placed. The constitu- 
tive parameters of effective fluid (p e for effective density and 
B e for effective bulk modulus) are unknown and to should 
be estimated. Pressure (^(r, (j>, t)) and radial component of 
velocity (v r (r, </>, t)) in host material and in effective material 
can be decomposed by cylindrical harmonics |9|: 
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where J n () and H n () are Bessel and Hankel functions of the 
first kind, n-th order, X™, Y™ are unknowns and m corre- 
sponds to e in effective media and for in host material. 

On the interface between the host fluid and the effective 
medium standard boundary conditions of velocity and pres- 
sure continuity should be satisfied. Due to cylindrical sym- 
metry and isotropy of materials, we could satisfy boundary 
conditions for each cylindrical harmonics independently |9): 
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where k e is the wave-number in effective medium, ko is the 
wave-number in host medium and e' m * and e~ 4Wt factors are 
dropped out. 

Main CPA condition is an absence of scattering on the 
boundary, which means that there is a perfect matching be- 
tween the coated inclusion and surrounding effective medium. 
Setting X ? c ; = 1 and Y£ = in (|2) we obtain the CPA equa- 
tion: 



X„ 



~ ^J' n (k Q L)J(k e L) + ^J n (k L)J'(k e L) = 
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(3) 



Exploiting the fact that metamaterial should be quasi-isotropic, 
in other words, the phase of the wave in host material should 
not change much on a distance between neighbor inclusions 
(koL —> and k e L —> 0). We use asymptotic formulas for 
Bessel and Hankel functions for near-zero argument limit. 
For small argument approximation only two, namely, zero-th 
and first orders of Bessel and Hankel functions are signifi- 
cant. After some algebra we obtain approximate formulas for 
effective constitutive parameters: 
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2.2 Surface acoustic wave 

Along the boundary between two ideal fluids any surface wave 
cannot propagate, because the boundary conditions cannot be 
satisfied. But if the one of materials has one negative con- 
stitutive parameter boundary states can be excited 1 1 1 ]. This 
is a direct analogy to surface plasmon states in plasma with 
£<0|. 

Let us consider a wave traveling on the interface between 
a half-space occupied by pure host fluid and a half-space oc- 
cupied by a fluid with embedded rods in it. Matching bound- 
ary impedances of evanescent waves in both composite and 
pure fluid, one can get: 



p e k z e - ppfcg 



= 0, 



(7) 



where fcg and k z are the projections of the wave- vector per- 
pendicular to the interface in a pure fluid and metamaterial, 
respectively. Dispersion equation for surface wave has the 
following form: 



(poPe){po/B e - p e /B ) 
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We show below, that the boundary states exist only if one of 
the constitutive parameters is negative. 



and D„ = Y n /X n is a scattering coefficient of the one em- 
bedded rod in an infinite fluid. In our case it is scattering 
of the acoustic wave in fluid on a round solid cylinder. This 
problem was widely investigated in past flO) . 

As it is seen from formulas Q, dynamical effective con- 
stitutive parameters are frequency dependent and in the case 
of resonance inside inclusion their values can be changed 
dramatically and even become negative. Thus the speed of 
acoustic wave inside the inclusions should be much smaller 
than in a host fluid, to form resonance frequency at which the 
wavelength in a host is much longer that distance between 
inclusions. 

Also we should notice, that bulk modulus depends on 
a monopolar scattering coefficient and effective density de- 
pends on the dipole coefficient. This is true, because the monopole 
mode is the uniform compression or expansion of the inclu- 
sion and dipole mode is the inclusion shift. 

The refractive index of the medium is described by the 
following formula: 



3 Results 



3.1 Bulk wave 

In a model used for calculations the rods are considered to 
be two-layered cylinders. Shell of the cylinder is soft sili- 
cone resirQ(> = 1.3 • 10 3 kg/m 3 , c n = 5.2 • 10 7 N/m 2 , 
c 44 = 3.25 • 10 6 N/m 2 , Loss: 600 dB/m at 1 MHz), core of 
the cylinder is poly-crystal silver JL?) (p = 10.5 • I0 3 kg/m 3 , 
cn = 1.397 • 10 n N/m 2 , c 44 = 2.7 • 10 10 A^/m 2 , Loss: 
40 dB/m at 1 MHz). Outer radius of the rod is 3 mm, ra- 
dius of the core is 0.5 mm, average distance between cylin- 
ders is taken 5 mm. Host material is water (p — 10 3 kg/m 3 , 
B = 2.25 • 10 9 iV/m 2 ). 

Fig. [T] shows results of calculations for bulk wave in the 
composite. At low frequencies effective constitutive parame- 
ters tend to be close to "classic" values for composites fPJ) . 
But at frequencies close to resonances inside the inclusions 
the values of dynamical effective constitutive parameters are 



(6) 



If the argument of refractive index tends to zero the wave- 
vector is pointing same direction as the Poynting vector. Ar- 
gument's value close to tt/2 means that the one of the consti- 
tutive parameters is negative and the other positive, thus the 
wave finally decays. Argument's value close to ir means that 
wave-vector and Poynting vector are directed oppositely and 
the wave have backward behavior. 



dramatically changed. In Fig. 1(a) two frequencies of monopo- 
lar resonances are seen and close to the resonances there are 
two bands where real part of compliance is negative (In this 
section we use compliance (s=\/B) instead of stiffness to 
obtain symmetry between constitutive parameters in |6])). 
In Fig. 1(b) frequency dependence of the dynamical ef- 



fective density is shown. Due to the complex shear and longi- 
tudinal field structure inside inclusion, there are two families 



1 There are various types of silicone resins, therefore our param- 
eters are approximate |4j 
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Fig. 1 Frequency dependence of effective compliance (a) den- 
sity |(b)| refractive index |(c)| and the argument of refractive index |(d)| 
for composite, normalized on zero-frequency limit 
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Fig. 2 Slowness (k s /ko) of surface acoustic wave at the interface 
between metamaterial and "normal" material 



of resonances associated with shear waves and longitudinal 
waves (S and L in the Figure, respectively). Thus dynamical 
effective density has a number of frequency ranges where real 
part is negative. The amplitude of resonances associated with 
shear components become greater because of inertia of the 
core. 

In Figs. 1 1(c)] and 1 1(d)] dependence of the refractive index 
and its argument are plotted. At the low frequency less than 
all resonance frequencies of inclusions composite exhibits 
properties of "normal" medium (region 1 ). At frequency greater 
that first monopolar resonance in the inclusions, the real part 
of 1/B become negative and the phase of wave vector is 
shifted to ir/2 (region 2). In this range wave-number is close 
to be pure imaginary, thus no oscillations can be exited and 
no wave is able to propagate. In such way a band gap in the 
spectrum is formed. On frequency greater than dipole reso- 
nance, effective density become negative and phase is shifted 
again for another ir/2 and the wave-number becomes close 
to be real negative. Thus phase velocity of the wave points 
backwards and opposite to the Poynting vector. Also, depen- 
dence of refractive index shows, that in the frequency band 
where constitutive parameters have simultaneous significant 
negative parts, the imaginary part of refractive index does not 
dominate and the wave is still able to propagate despite of 
losses in inclusions. 



3.2 Surface acoustic wave 

Fig.|2]shows slowness (k e /ko) for surface wave to bulk wave 
in water. As it was predicted, wave-number does not have 
real part at low frequencies. Since the parameters of meta- 
material are positive, the wave cannot propagate along the 
interface. But at frequency greater than resonance (compare 



to Fig. 1(a) i one of material parameters become negative and 
wave-number is real and the surface states appear. Surface 
wave is much slower than bulk wave, thus the surface states 
are bounded to the interface. 



3.3 Complex Doppler effect 



It was shown that in electromagnetic media with simultane- 
ous negative e and \i Doppler effect should be reversed (T]. 
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Fig. 3 Graphical solution for equation {9} in metamaterial 

Similar effect can exist for acoustic metamaterials. Let us 
consider a source travels inside a composite with velocity V. 
It transmits signal with an angular frequency ui s . We can as- 
sume that source is moving in one line with receiver. To find 
frequency on the receiver (w) we should solve the following 
equation: 

uo s = u - V ■ k(uj) (9) 

For dispersive medium like metamaterial solution for equa- 
tion Q is not trivial fl4) . 

Graphical solution for equation |9| is plotted in Fig. [3] 
We assume that dumping is small and and wave-number is 
real. Angular frequency of source is taken 1.22 • 10 5 rad/s 
and speed is 15 m/s. The source moves away from receiver. 
For the pure host material the result is predictable, the fre- 
quency on receiver is down-shifted. In contrast, for compos- 
ite equation |9| has two roots and the two modes with down- 
shifted and up-shifted frequencies are present. Generation of 
two modes is possible, because modes have different phase 
velocities and one of the waves if forward and another is 
backward (14). 
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4 Conclusions 

In this paper effective dynamic density and stiffness of com- 
posite consisted of water with embedded solid two-layered 
rods in terms of coherent potential approximation are cal- 
culated. It is shown, that there are some frequency ranges 
in which dynamic constitutive parameters are simultaneously 
negative and the propagating wave becomes backward. 

Dispersion of the surface acoustic wave at the interface 
between metamaterial and "normal" material is calculated. 
It is shown that there are frequency ranges in which surface 
states are bounded to the interface. 

Doppler effect in composite is considered. It is shown, 
that in metamaterial Doppler shift is reversed and several modes 
could be excited. 
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